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PREFACE 


The  investigation  reported  herein  was  conducted  for  the  Office, 
Chief  of  Engineers,  U.  S.  Army,  by  personnel  of  the  Soil  Dynamics  Divi- 
sion (SDD),  Soils  and  Pavements  Laboratory  (S&PL),  U.  S.  Army  Engineer 
Waterways  Experiment  Station  (WES),  as  a part  of  Project  HAl6ll02AT?2 , 
Task  A2,  Work  Unit/06,  "Effectiveness  of  Earth  Penetrators  in  Various 
Geologic  Environments." 

Mr.  R.  S.  Bernard  conducted  the  research  during  the  period 
October  1977-January  1978  with  the  technical  guidance  of  Dr.  B.  Rohani 
€nl  under  the  supervision  of  Dr.  J.  G.  Jackson,  Jr.,  Chief,  SDD,  and 
Messrs.  J.  P.  Sale  and  R.  G.  Ahlvin,  Chief  and  Assistant  Chief,  S&PL, 
respectively.  This  report  was  also  prepared  by  Mr.  Bernard. 

COL  J.  L.  Cannon,  CE,  was  Director  of  WES  throughout  the  investi- 
gation and  during  the  preparation  of  the  report.  Mr.  F.  R.  Brown  was 
Technical  Director. 
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DEPTH  AMD  MOTION  PREDICTION  FOR  EARTH  PENETRATORS 


PART  I:  INTRODUCTION 


1.  Earth-penetrating  weapons  (EPW's)  offer  a means  of  reducing 
collateral  effects  in  the  selective  destruction  of  localized  targets 
(airfields,  factories,  utilities,  etc.).  The  effectiveness  of  these 
weapons,  however,  is  contingent  upon  (a)  accurate  delivery,  (b)  impact 
survival,  and  (c)  proper  depth-of-burst  (DOB).  The  first  two  require- 
ments are  obvious;  the  last  is  important  because  there  is  an  optimum 
DOB  in  the  trade-^ff  between  destructive  capability  and  collateral 
effects.  To  guarantee  detonation  at  or  near  the  optimum  DOB,  it  is 
necessary  to  be  able  to  predict  the  subsurface  motion  of  an  EPW. 

2.  There  are  three  alternatives  for  predicting  EFW-motion  after 
impact;  (a)  two-dimensional  finite-difference  codes,  (b)  equations  of 
motion  based  partly  on  theory,  and  (c)  empirical  formulae.  Alternative 
(a)  is  acceptable  in  some  cases  but  is  too  expensive  for  general  appli- 
cation. Alternatives  (b)  and  (c)  require  far  less  computation  and  are 
therefore  suitable  for  parameter  studies  and  multiple  predictions. 

3.  Among  the  theoretically  based  equations  of  motion,  those 
developed  from  the  Cavity  Expansion  Theory  (CET)  have  been  preferred 
at  the  U.  S.  Army  Engineer  Waterways  Experiment  Station.  Of  the  purely 
empirical  formulae.  Young's  equation**-^  seems  now  to  be  the  most  widely 
accepted  for  calculating  final  penetration  depth  in  soil.  The  main  dif- 
ference between  these  two  approaches  lies  in  the  target  description.  The 
CET-based  analysis  calculates  penetration  resistance  in  terms  of  the  me- 
chanical properties  of  the  target  (density,  strength,  elasticity,  and 
compressibility).  Young's  equation,  on  the  other  hand,  employs  a single 
empirical  parameter  (S)  to  quantify  the  target  penetrability: 
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Table  1* 

Typical  S-Numbers  for  Natural  Earth  Materials 


0.2-1 


1-2 


2-3 


h-6 


8-12 

10-15 

20-30 

1*0-50 


Materials 


Massive  medium-  to  high-strength  rock,  with  few  fractures. 
Concrete,  2000  to  5000  psi,  reinforced. 

Frozen  silt  or  clay,  saturated,  very  hard.  Rock,  weathered, 
low  strength,  fractured.  Sea  or  freshwater  ice  more  than 
10  feet  thick. 

Massive  gypsite  deposits  (WSMR**)-  WeJJL-cemented  coarse  sand 
and  gravel.  Caliche,  dry.  Frozen  moist  silt  or  clay. 

Sea  or  freshwater  ice  from  1 to  3 feet  thick.  Medium  dense, 
medium  to  coarse  sand,  no  cementation,  wet  or  dry.  Hard,  dry 
dense  silt  or  clay  (TTRt  dry  lake  play as ):„ ! Desert  alluvium. 

Very  loose  fine  sand,  excluding  topsoil.  Moist  stiff  clay  or 
silt,  medium  dense,  less  than  about  50  percent  sand. 

Moist  topsoil,  loose,  with  some  clay  or  silt.  Moist  medium 
stiff  clay,  medium  dense,  with  some  sand. 

Loose  moist  topsoil  with  humus  material,  mostly  sand  and  silt. 
Moist  to  wet  clay,  soft,  low  shear  strength. 

Very  loose  dry  sandy  topsoil  (Eglin  AFB).  Saturated  very 
soft  clay  and  silts,  with  very  low  shear  strengths  and  high 
plasticity.  (Great  Salt  Lake  Desert  and  bay  mud  at  Skaggs 
Island.)  Wet  lateritic  clays.  j 


* Taken  from  Reference  6. 

**  White  Sands  Missile  Range,  New  Mexico, 
t Tcnopah  Test  Range,  Nevada. 


Table  2* 

Rose  Performance  Coefficient 


Hose  Shrpe 

Nose 

Length-to-Diameter 
Ratio  (L/D ) 

N 

Flat 

0 

0.56 

Hemisphere 

0.5 

0.65 

Cone 

1 

0.82 

Tangent  ogive** 

l.U 

0.82 

Tangent  ogive 

2 

0.92 

Tangent  ogive 

2.k 

1.0 

Inverse  ogive 

2 

1.03 

Cone 

2 

1.08 

Tangent  ogive 

3 

1.11 

Tangent  ogive 

3.5 

1.19 

Step  cone 

3 

1.28 

Biconic 

3 

1.31 

Inverse  ogive 

3 

1.32 

Cone 

3 

1.33 

* Taken  from  Reference  6. 

**  For  tangent  ogives,  L/D  is  related 
to  the  caliber  radius  (CRH)  by 
CRH  = L2/D2  + 1/U  . 
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H.  A recent  study  by  Rohani  et  al.  has  shown  that,  for  in- 
accessible targets,  analysts  with  different  backgrounds  are  more  likely 
to  agree  upon  S-number  estimates  than  upon  mechanical-property  estimates. 
Equating  analyst  agreement  with  reduced  uncertainty  of  prediction. 

Young's  equation  is  apparently  better  suited  for  inaccessible  targets 
than  is  the  CET  analysis.  This  value  judgement  does  not  necessarily 
apply  for  accessible  targets,  where  S-numbers  and  mechanical  properties 
can  be  determined  by  experiment. 

5.  Young's  equation  has  been  fairly  well  v&lidated  for  soil,  but 

its  applicability  is  questionable  for  hard  materials,  such  as  concrete 

8 

and  rock.  In  fact,  penetration  data  for  concrete0  indicate  that  the 
final  depth  is  more  nearly  proportional  to  W/A  than  to  /W/A  . This 
seems  to  warrant  a composite  prediction  technique,  with  distinct 
analyses  for  soil  and  rock  (or  rocklike  materials). 

Purpose 

6.  An  improved  method  of  earth-penetration  -nalysis  is  sought, 
particularly  for  inaccessible  targets.  A projectile  equation  of  motion 
will  be  formulated  for  soil,  retaining  Young's  £>-number  and  approx- 
imately reproducing  Young's  final-depth  equation  after  integration.  A 
distinct  equation  of  motion  will  be  developed  for  rock,  using  density, 
unconfined  compressive  strength,  and  Rock  Quality  Designation  (RQD)^ 

to  describe  the  target.  The  two  analyses  will  then  be  interfaced  in 
a composite  analysis  for  targets  containing  both  soil  and  rock  layers. 

Scope 

7.  Part  II  contains  the  soil  penetration  analysis,  and  Part  III 
the  rock  penetration  analysis.  Equations  for  layered  targets  are  de- 
veloped in  Part  IV.  Part  V states  the  conclusions  drawn  from  the 
investigation. 
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PAST  II:  SOIL  PENETRATION  ANALfSIS 


8.  In  soil  penetration  tests  conducted  with  instrumented  projec- 
tiles, the  deceleration  curve-shape  most  often  observed  is  approximately 
that  of  a "square  puloe."^  In  targets  with  distinct  soil  layers,  the 
deceleration  record  is  usually  a series  of  twc  or  more  square  pulses, 
as  shown  in  Figure  2 (taken  from  Reference  10).*  The  magnitude  of  the 
deceleration  is  velocity  dependent;^  an  increase  in  the  impact  velocity 
usually  causes  a proportionate  increase  in  the  deceleration  level. 


DECELERATION.  «*• 

SOIL  PROFILE  0 100  200  SOO  400  SOO 


* The  oscillations  in  the  record  are  thought  to  arise  from  the  projec- 
tile nonrigidity  and  from  the  accelerometer  mountings  themselves. 
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9.  In  order  to  be  generally  consistent  vita  experimental  data, 
the  projectile  equation  of  motion  must  (a)  exhibit  velocity  dependence 
and  (b)  generate  flat  deceleration  records*  in  uniform  targets.  The 
simplest  equation  of  motion  that  can  satisfy  both  criteria  is 

-m  ft  ■ * m § * bv  + cz  (3) 

vhere 

M * projectile  mass 
v * instantaneous  velocity 
t * time 

z * instantaneous  depth 

b,c  * coefficients  dependent  or.  soil  type  and  projectile  param- 
eters but  independent  of  velocity  and  depth 

The  initial  conditions  for  Equation  3 are  v * V and  z * 0 ; the  final 
conditions  are  v * 0 and  z * Z . 

10.  Disregarding  any  change  in  load  during  the  nose-embedment 
process.  Equation  3 is  integrated  with  respect  to  z t obtaining  the 
following  expression: 

« cZ2  ♦ 2b 

It  can  be  shown  that  if  dv/dt  Is  approximately  constant  between  z * 0 
and  z * Z , then  Equation  U reduces  to 

MV2  * cZ2  + | bVZ  (5) 

Thus,  the  expression  for  the  final  penetration  depth  becomes 

2 “ c ("  f b *2  + McJ  (6) 


* "Flat"  in  the  sense  tnat  the  deceleration  is  essentially  constant 
after  nose  embedment. 
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Since  the  deceleration  curve  is  nearly  flat,  the  initial  deceleration 
(at  z =0)  is  about  equal  to  the  final  deceleration  (at  v = 0),  i.e., 


bV  = cZ 


(7) 


Combining  Equations  6 and  7 to  eliminate  V and  Z , it  follows  that 

b « - |b  + yj^b2  4 Me  (8) 

Upon  rearranging  and  squaring  both  sides  of  Equation  8,  a simple  rela- 
tion occurs  between  b and  c : 

b2  = y Me  (9) 

Equation  9 reveals  that  b and  c are  not  independent;  b is  pro- 
portional to  the  square  root  of  the  projectile  mass: 

b = ^Tmc  (10) 

Thus,  the  presumption  of  a flat  deceleration  curve  requires  that  b 
depend  on  M , while  c remains  independent  of  M . * 

11.  Sandia  Laboratories  has  conducted  at  least  bj  low-speed 
penetration  tests  (V  < 80  m/sec)  in  the  TTR  Main  ^ake  area  . In  *3 
of  these  tests,  the  projectiles  had  nearly  the  same  mass  (i00  to 
105  kg),  with  different  diameters  and  nose  shapes.  The  trends  in  the 
TTR  data  (and  likewise  in  Young's  equation)  can  be  summarized  as 
follows : 

£i.  The  relation  between  final  depth  (Z)  and  impact  velocity 
(V)  is  nonlinear  when  V is  significantly  less  than 
60  m/sec. 


* It  might  be  argued  that  M should  appear  only  on  the  left-hand  side 
of  the  projectile  equation  of  motion  (Equation  3).  The  presumption  of 
a flat  deceleration  curve  at  the  outset,  however,  forces  M to  appear 
on  the  right-hand  side  in  the  velocity  coefficient  b . No  additional 
interpretation  can  be  given  to  the  mass  dependence  of  b without  ex- 
plaining why  the  deceleraticn  curve  is  flat  in  the  first  place. 
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t_.  The  final  depth  is  inversely  proportional  to  the  prcjec- 
. tile  diameter  (D). 

£.  The  final  depth  is  a function  of  the  projectile  rose 
shape. 

The  ncnlinear  relation  between  Z and  V at  low  speed  suggests  that 
the  equation  of  notion  needs  an  additional  constant  tern: 


u dv  ..  dv  _ , 

-K  — — = -Mv  — * a + bv  + cz 
dt  az 


(11, 


With  this  modification  the  deceleration  curve  still  remains  flat,  and 
Equation  11  can  be  solved  for  the  final  depth  in  the  same  way  as 
Equation  3.  The  solution  obtained  is 


7 t i 


m 

« 

-(a  + | bv)  + 

| a + f bvj2  * XcV2 

( ' o > 


Assuming  that  Z is  inversely  proportional  to  D and  directly  prcpcr- 
tional  to  X and  £ , it  fellows  from  Equations  1C  and  12  that 


a; 


b * 


* 'L\ 


e:‘ 


£“K£ 


where  a and  £ are  universal  constants,  independent  of  projectile 


parameters  and  so.il  properties. 

12.  Up  to  this  point,  the  analysis  has  dealt  only  with  generali- 
ties, establishing  the  form,  cf  the  equation  of  motion  (Equation  11/  and 
the  resulting  form  of  the  final-depth  equation  (Equation  12}.  The  r.cse- 
perfornance  coefficient  (K)  and  the  penetrability  index  (S)  have  been 
retained  for  consistency  with  Young's  equation.  The  numerical  values 
cf  u ar.d  £ must  now  be  obtained  by  fitting  Equation  12  to  actual 
penetration  data.  The  best  fit  to  the  TTE  data  (Figure  3)  occurs  for 
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Figure  3.  Correlation  of  TTR  penetration  data 
with  Equation  12  for  empirical  determination 


of  a 

and 

e 

a * 2.2 

x io6 

N/m 

(16) 

6 = 2.8 

x 107 

K/m3 

(17) 

where  Young’s  value  of  S * 5.2  has  been  assumed  for  the  TTR  Main 
Lake  area.^  With  a and  6 given  in  the  units  shown,  Z can  be  ob- 
tained in  metres  from  Equation  12  (in  which  the  values  of  M , D , 
and  V must  be  expressed  in  kilograms,  metres,  and  metres  per  second, 
respectively). 

13.  The  soil  penetration  analysis  is  now  complete.  Due  to  the 
presumption  of  a flat  deceleration  curve  in  the  formul;icion  of  the 
equation  of  motion,  the  effect  of  the  projectile  mass  (or  weight)  in 
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Equation  12  is  slightly  different  from  the  effect  in  Young’s  equation. 
Close  inspection  of  Equation  12  reveals  that  Z * M as  MV^  •*>  0 and 
and  that  Z « as  MV^  -*•  » . As  the  impact  energy  increases  from 
zero,  the  M-dependence  of  Z gradually  changes  from  linear  to  square 
root.  In  Young's  equation,  the  M-dependence  changer  from  linear 
(at  M a,  l kg)  to  square  root  (at  M >_  27  kg),  indepen.  nt  of  the 
impact  velocity.  Figure  1*  shows  a comparison  of  Equation  12  with 
Young's  equation  (for  W ^ 27  kg)  and  with  normalized  soil  penetration 
data. 

ll*.  Figure  5 shows  a comparison  of  calculated  deceleration 


Figure  U.  Comparison  of  Equation  12  with 
Young's  equation  and  with  normalized  soil 
penetration  data 
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Figure  5-  Comparison  of  calculated 
deceleration  records  in  hypothetical 
soil  target  for  impact  velocity  of 
100  m/sec 


curves*  for  a 500-kg,  25-cm  projectile  and  a 5-kg,  2.5-cm  projectile  at 
an  impact  velocity  of  100  m/sec.  Figure  6 shows  a similar  comparison 
at  300  m/sec.  In  both  cases,  the  lighter  projectile  incurs  the  higher 
deceleration,  but  the  difference  in  deceleration  between  the  two  pro- 
jectiles decreases  as  V increases. 


* These  results  were  obtained  by  numerical  integration  of  Equation  11 
accounting  for  the  change  in  diameter  during  nose  embedment.  The 

projectiles  have  the  same  value  of  W/A  . 


PART  III:  ROCK  PENETRATION  ANALYSIS 


15.  Reference  12  documents  an  empirical  equation  for  calculating 
final  penetration  depth  in  massive  rock  deposits.  Expressed  in  dimen- 
sionless form,*  the  equation  is: 


eL. 

M/A 


0.2V 


/r/joo\°-8 
V * \R  QD/ 


(18) 


where 

p * mass  density  of  the  rock 

Y <=  unconfined  compressive  strength  of  the  intact  rock 

9 

The  RQD,  originally  proposed  by  Deere,  is  an  index  for  the  degree  of 
fracturing  of  the  rock  (in  situ)  at  a given  site.  Its  value  (in  per-  ' 
cent)  is  determined  hy  a special  core-logging  procedure: 

All  solid  pieces  of  core  that  are  10  cm  long  or  longer 
are  added  up,  and  this  length  is  called  the  modified 
core  recovery.  The  modified  core  recovery  is  divided 
by  the  total  length  of  core  run,  and  the  quotient  mul- 
tiplied by  100  percent  is  the  value  of  the  RQD. 

16.  Equation  18,  like  Young's  equation,  gives  only  the  final 

depth  of  penetration;  it  does  not  predict  the  projectile  motion  after 

impact.  This  equation  was  obtained  by  curve- fitting  data  from  rock 
12 

penetration  tests  in  which  p , Y , and  RQD  were  all  known.  The 
data  were  somewhat  scattered,  however,  and  the  linearity  in  V was 
chosen  mainly  for  convenience  of  calculation. **  A linear  equation 
seemed  to  fit  the  overall  data  at  least  as  accurately  as  other  func- 
tions of  V . 


* Obviously,  the  values  of  p,  Z , M,  A,  Y,  and  V must  be  ex- 
pressed in  compatible  units  if  Equation  18  is  to  be  dimensionless. 

**  Most  of  the  scatter  is  due  to  uncertainty  in  the  coefficient  of 
penetrability  from  one  target  to  the  next.  When  several  data  exist 
for  a single  target,  they  usually  exhibit  (a)  less  scatter  and  (b)  a 
weakly  nonlinear  relation  between  Z and  V . However,  when  data 
from  multiple  targets  are  superposed,  the  apparent  nonlinearity  in  V 
is  lost  in  the  scatter. 
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17.  Existing  deceleration  record  are  too  few  to  allow  definitive 

statements  about  the  shape  of  the  deceleration  curve  in  rock.  Thus, 

the  assumption  of  a flat  deceleration  record  is  not  Justified  for  rock 

as  it  was  for  soil.  The  formulation  of  the  projectile  equation  of 

motion  must  rest  on  final-depth  data  alone. 

8 

18.  Canfield  and  Clator  have  conducted  a series  of  penetration 
tests  in  concrete,  which  is  similar  in  some  ways  to  intact  (unjointed, 
unfractured)  rock.  The  test  results,  obtained  for  7 .6-  and  76-mm 
projectiles,  suggest  the  following  conclusions: 

a.  The  final  depth  is  directly  proportional  to  the  pro- 
jectile mass  and  inversely  proportional  to  the  square 
of  the  projectile  diameter. 

b.  The  final  depth  increases  linearly  with  impact  velocity 
for  V > 300  m/sec,  but  the  curve  has  a nonzero  in- 
tercept, indicating  nonlinearity  at  low  velocity. 

A two-term  equation  of  motion  is  needed  to  reproduce  the  trends  in  the 
concrete  penetration  data: 

-M  ||  * -Mv  g * J D2(a'  + b’v)  (19) 

where  a'  and  b'  are  coefficients  independent  of  v and  z but 
dependent  on  target  properties.  Neglecting  the  change  in  D during 
the  nose— embedding  process,  the  solution  for  the  final  depth  is 


Presuming  that  rock  and  concrete  respond  in  much  the  same  way  to 
penetration.  Equations  19  and  20  should  also  be  applicable  for  rock, 
though  the  values  of  a*  and  b'  may  be  different  from  those  for 
concrete . 

19.  Figure  7 shows  the  concrete  penetration  data  obtained  by 
Canfield  and  Clator  for  a 76-mm  projectile.  Equation  22  fits  these 
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Figure  7.  Correlation  of  penetration  data  for 
5000-psi  concrete  vith  Equation  20  for  empirical 
determination  of  a'  and  b' 


data  fairly  well  when  a'  and  b*  are  given  the  following  values:* 

(a*  « 1.6Y  (21) 

concrete:  < 

: |,t  ,3.6,/pY  - - (22) 

In  the  Canfield  and  Clator  tests,  the  unconfined  compressive  strength 
8 o o 

was  3.1*5  x 10  dyne/cm  , and  the  density  2.88  gm/cm  . 

20.  Equations  19  and  20  are  applicable  for  rock  penetration  (to 

at  least  the  same  degree  as  Equation  18)  when  Equations  21  and  22  are 

amended  as  follows: 


* The  values  jhosen  for  a’  and  . * do  not  produce  a very  good  fit 
to  the  concrete  data  as  such,  buc  they  do  produce  a good  fit  to 
superposed  concrete  and  rock  data. 
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Figure  8 shows  a (nondimensional)  comparison  of  Equations  18  and  20  with 
the  rock  penetration  data  used  in  the  formulation  of  Equation  18. 


SYMBOL 

TARGET 

ROD 

average 

strength 

BARS 

AVERAGE 

DENSITY 

3m/c*> 

A 

WCLOEO  TUEE 

100 

•00 

1.BB 

a 

SANOSTONE 

•2 

234 

2.00 

A 

WELDED  AGGLOMERATE 

*0 

275 

1.02 

• 

SANOSTONE 

S7 

«** 

2.12 

■ 

SANDSTONE 

33 

4CO 

2.(4 

f 

GRANITE 

U 

M2 

2.02 

penetration  equations  with  rock  penetration 
data 


21.  Figure  9 shows  calculated  deceleration  curves*  in  concrete 


* These  curves  were  obtained  by  numerical  integration  of  Equation  19, 
allowing  for  the  change  in  diameter  during  nose  embedment. 
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Figure  9«  Calculated  deceleration  curves  for 
5000-psi  concrete 


for  three  different  impact  velocities.  Due  to  the  form  of  Equation  19 
and  the  relative  magnitudes  of  a’  and  V , the  maximum  deceleration 
occurs  at  the  point  of  nose  embedment,  vith  the  minimum  deceleration  at 
the  final  penetration  depth.  If  the  coefficient  of  Y (Equation  21) 
vere  larger,  and  the  coefficient  of  (Equation  22)  sma?ler,  then 

the  calculated  deceleration  curves  vould  be  flatter,  and  the  relation 
between  Z and  V more  nonlinear. 
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PART  IV:  LAYERED  TARGET  ANALYSIS 


22.  Few  natural  targets  are  uniform  to  any  great  depth,  and  it  is 
not  uncommon  for  earth  penetrators  to  encounter  abrupt  changes  in  pene- 
trability. This  poses  no  special  problem  if  the  adjacent  layers  are 
similar  materials  (both  r.oil  or  both  rock),  but  a difficulty  arises 
when  a layer  of  rock  overlies  a leyer  of  soil.  A finite  thickness  of 
rock  is  more  penetrable  (per  unit  depth)  than  a half-space  of  the  same 
rock,  yet  there  do  not  exist  enough  rock-over-soil  penetration  data  to 
formulate  a quantitative  model  for  tne  weakening  of  the  rock  near  the 
interface.  Aside  from  this  shortcoming,  the  penetration  analysis  for 
layered  targets  is  a logical  extension  of  the  work  in  Parts  II  and  III. 


Equation  of  Motion  for  Layered  Targets 


23.  Upon  considering  the  projectile  equation  of  motion  (Equa- 
tions ll  and  19),  it  is  convenient  to  introduce  the  quantity  o , which 
represents  the  resisting  force  per  unit  frontal  area.  This  allows 
Equations  11  and  19  to  be  rewritten  as 


,,  dv  k _2„ 
”M  dt  “ TT  D ° 


(25) 


Whenever  the  projectile  nose  is  in 'contact  with  two  adjacent  layers 
(Figure  10),  Equation  25  is  replaced  by 


’2ff2 


(26) 


where 


D Dp  — maximum  diameter  in  contact  with  layers  1 and  2, 
respectively 

c^,Og  = ,ralue  of  o in  layers  1 and  2,  respectively 
1 Layered  Soil  Targets 

2U.  Equation  26  is  applicable  for  soil  and  rock  layers  alike,  so 
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Figure  10.  Projectile  in  contact  with 
two  layers  at  the  same  time 


long  as  the  definitions  of  and  make  Equation  26  consistent  with 

Equations  11  and  19.  For  soil  layers  in  particular,  and  must 

be  defined  such  that  Equation  26  reduces  to 

"M  dt  = U1  + a2)  + (bl  + VV  + (C1  + C2)Z  (27) 


where 


1 


I 


♦ 


and  and  are  the  S-numbers  in  layers  1 and  2,  respectively. 

With  the  coefficients  defined  in  Equations  28-33*  Equation  2T  collapses 
to  Equation  11  whenever  ^ In  order  for  Equation  26  ’ o be 

consistent  (for  soil)  with  Equation  2T,  it  then  follows  that 


where  z is  the  depth  of  the  nose  tip.  r*asured  from  the  target  surface. 
Figure  11  shows  a numerically  calculated  deceleration  curve  for  a hypo- 
thetical three-layer  soil  target. 


Layered  Rock  Targets 


25.  In  order  to  make  Equation  26  consistent  (for  rock)  with 
Equation  19*  the  definitions  required  for  and  o g are: 
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crh  - e.o 
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V=  100  r/i«c 


Figure  11.  Calculated  deceleration  record  for  a 
hypothetical  three-layer  soil  target 

vhere  the  subscripts  1 and  2 designate  the  evaluation  of  p , Y , and 
RQD  in  layers  1 and  2,  respectively.  Figure  12  shows  a numerically 
calculated  deceleration  record  for  a hypothetical  three-layer  rock 

target . 

Composite  Targets 


26.  Equation  26  is  applicable  for  targets  containing  both  soil 
and  rock  layers.  Equations  3k-35  and  36-37  define  o in  the  soil  and 
rock  layers,  respectively.  The  tc-tal  force  on  the  projectile  can  be 
obtained  by  substituting  the  appropriate  expressions  for  and 

into  Equation  26.  For  example,  if  layer  1 is  soil  and  layer  2 is  rock, 
then  and  are  given  by  Equations  3k  and  37,  respectively. 

Figure  13  shows  a numerically  calculated  deceleration  curve  for  a hypo- 
thetical soil/rock/soil  target. 
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Figure  12.  Calculated  deceleration  record  for  a 
hypothetical  three-layer  rock  target 


Figure  13.  Calculated  deceleration  record  for  a 
hypothetical  soil/rock/soil  target 


PART  V:  CONCLUSIONS 


27.  The  soil  penetration  analysis  (Part  II)  is  a reexamination  of 
the  data  originally  used  in  the  formulation  of  Young's  equation.  The 
projectile  equation  of  motion  (Equation  11)  obtained  therefrom  (a)  offers 
a means  of  calculating  the  projectile  motion  explicitly  and  (b)  appears 
to  clarify  the  effect  of  the  projectile  mass  upon  the  final  penetration 
depth.  The  accuracy  of  prediction  is  the  same  as  for  Young's  equation: 
rigid-body  deceleration  and  final  depth  +20  percept,  when  the  S-number 

is  known  accurately. 

28.  The  rock  penetration  analysis  (Part  III)  is  based  on  fewer 
jdata  than  the  3oil  penetration  analysis.  Nevertheless,  for  high-quality 
| rock  (RQD  > 90)  in  which  the  strength  is  known  accurately,  the  final- 

| depth  predictions  are  usually  accurate  within  +20  percent.  For  the 
same  situation,  predictions  of  peak  rigid-body  deceleration  should  be 

I accurate  at  least  within  +50  percent.  In  order  to  improve  upon  this, 

) 

| more  projectile  deceleration  records  are  needed  to  get  a better  under- 
I standing  of  (a)  the  general  shape  of  the  deceleration  curve  and  (b)  the 

! 

[relation  between  rock  properties  and  rock  penetrability. 

29.  The  layered  target  analysis  (Part  IV)  is  a logical  extension 
[of  Parts  II  and  III.  The  equations  therein  should  be  about  as  accurate 

| for  layered  targets  as  for  unlayered  targets , with  the  exception  of  rock 
layers  over  soil.  In  rock  layers,  the  resistance  to  penetraticn  de- 
creases near  a rock/soil  interface,  but  the  analysis  does  not  account 
for  this  directly.  For  calculations  involving  thick  rock  layers  over 
soil,  the  reduced  penetrability  can  be  approximated  indirectly,  however, 
by  using  a simple  rule  of  thumb:  Reduce  the  rock  thickness  by  two  or 
three  projectile  diameters,  and  increase  the  thickness  of  the  under- 
lying soil  by  the  same  amount. 
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APPENDIX  A:  NOTATION 


a,b,c  Coefficients  in  projectile  equation  of  motion  for  soil 
(Equation  11 ) 

a'  ,b'  Coefficients  in  projectile  equation  of  motion  for  rock 
(Equation  19) 

2 

A Projectional  cross-sectional  area  (ttD  A) 

CRH  Ogive  caliber  radius 

D Projectile  diameter 

g Gravitational  acceleration  (9.8  m/sec^") 

K Young's  mass-scaling  factor 

L Projectile  nose  length 

M Projectile  mass 

N Projectile  nose-performance  coefficient 

RQD  Rock  Quality  Designation 

S Young's  penetrability  index  for  soil 

t Time 

v Instantaneous  velocity 

V Impact  velocity,  m/sec 

W Projectile  weight,  kg 

Y Unconfined  compressive  strength  for  rock  and  concrete 

z Instantaneous  depth 

Z Final  penetration  depth,  m 

a, 8 Coefficients  in  projectile  equation  of  motion  for  soil 

(Equations  11,  13,  lU,  15) 

ir  3.A16 

p Mass  density  for  rock  end  concrete 
o Resisting  force  per  unit  frontal  area 


Note:  Subscripts  1 and  2 denote  evaluation  of  a given  quantity  in 
layers  1 and  2,  respectively  (Figure  10 ). 
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